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Abstract. A two-dimensional bistable lattice is a periodic triangular network of non-linear bi- 
stable rods. The energy of each rod is piecewise quadratic and has two minima. Consequently, a 
rod undergoes a reversible phase transition when its elongation reaches a critical value. We study an 
energy minimization problem for such lattices. The objective is to characterize the effective energy 
of the system when the number of nodes in the network approaches infinity. The most important 
feature of the effective energy is its "fiat bottom". This means that the effective energy density is 
zero for all strains inside a certain three-dimensional set in the strain space. The flat bottom occurs 
because the microscopic discrete model has a large number of deformed states that carry no forces. 
We call such deformations still states. In the paper, we present a complete characterization of the 
"fiat bottom" set in terms of the parameters of the network. This is done by constructing a family 
of still states whose average strains densely fill the set in question. 

The two-dimensional case is more difficult than the previously studied case of one-dimensional 
chains, because the elongations in two-dimensional networks must satisfy certain compatibility con- 
ditions that do not arise in the one-dimensional case. We derive these conditions for small and 
arbitrary deformations. For small deformations a complete analysis is provided. 

Key words. Phase transformation in solids, Bi-stable lattice, Eigenstrain, Compatibility con- 
ditions. 
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1. Introduction. This paper suggests a simple discrete model for phase transi- 
tion theory. In two dimensions, we study deformation of bi-stable lattices (BL), which 
are periodic triangular networks of bistable elastic rods. A similar network with lin- 
early elastic rods has been used by Cauchy, who formally averaged this assembly and 
obtained one of the first equations of continuum elasticity [8] . Here, we consider the 
same assembly, assuming, however, that the rods are bi-stable. The bi-stability mod- 
els phase transitions: it is assumed that each rod has two equilibrium states that of 
length I and l{l + s), respectively. Here, s can vary between ^ and 2, so that a triangle 
of transformed rods still remains a triangle. 

We believe that this bi-stable lattice model is the simplest and natural finite di- 
mensional model for phase transitions in solids. Although it does not have all the 
features of continuum model, it appears to capture most of the expected features. At 
the same time this model allows for detailed description of inhomogeneous deforma- 
tions and does not involve ad hoc assumptions of the continuum theory. 

One dimensional chains of bi-stable rods are well investigated [TUl [TTl IT^ ITB] . 
Each bistable element can be in one of the two states {short and long mode) that differ 
by the equilibrium length. It is assumed for definiteness that initially all rods are in 
short mode and can transit to the long mode. The state of a chain with some transited 
elements is characterized by a scalar relative elongation, or one-dimensional strain, 
that is the normalized difference between the length in the deformed configuration 
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and reference length. 

Some results on the two-dimensional bi-stable lattices can be found in [T^l fT5] 
(explicit construction of certain special solutions, discrete Green functions), and [HIIT] 
(direct numerical simulations of time-dependent impact fracture problems). Other 
related problems are treated in the book 14J that also contains a wealth of additional 
references. 

The present work differs from the above mentioned papers both in approach 
and in scope. It addresses the specific difficultly of the two-dimensional problem 
- compatibility conditions on the rods' length. These conditions are automatically 
satisfied in one dimensional chains. In two-dimensional lattice, the transition of links 
to a new state generally leads to elastic deformation of all rods. For example, one 
transformed (elongated) rod does not fit into a triangular lattice of non-transformed 
(short) links without distorting the lattice. Each such distortion increases the energy 
of the assembly. 

One of the objectives of the paper is to describe special lattice deformations that 
we call still states. These deformations carry zero forces in each rod, and thus the 
energy of each still state is zero. In any such state the length of each rod should be 
either I or Z(l -I- s). Since the number of rods is larger than the number of nodes times 
the dimension, an arbitrary assignment of lengths may not correspond to an actual 
lattice deformation. Therefore compatibility conditions must be satisfied at a still 
state. 

We describe the set of average strains of still states. First, we consider the network 
of rods undergoing arbitrary deformations. We assume that an effective strain tensor 
is represented by its eigenvalues. In the plane of eigenvalues, the set Hl of all still 
states is a nonconvex curvilinear hexagon. We demonstrate structures that densely 
fill Ti.L and derive nonlinear compatibility conditions. 

Further, we consider the case of small deformation, linearize the compatibility 
conditions and develop the theory in more detail. We show that the linearized com- 
patibility conditions are necessary and sufhcient for existence of a lattice deformation 
with a given set of link elongations. In this sense, the conditions are analogous to the 
compatibility conditions in continuum elasticity. We also describe a sequence of still 
states whose average strains reach the average strain of any still state by incremental 
change of the structure starting from the reference configuration. This procedure can 
describe an evolutionary irreversible process of phase transition, such as the damage 
distribution. In this case, the initial state corresponds to "undamaged" configuration 
and structural change - to irreversible damage propagation. 

We also construct the set V of all average strain tensors corresponding to still 
states and sets of deformed elements that densely fill that set. In the space of entries 
of the strain tensor, the set V turns out to be a parallelepiped. We obtain explicit 
formulas that relate the side lengths and orientation of this parallelepiped to the 
non-dimensional critical elongation parameter s, and three lattice direction vectors. 

Finally, we suggest an equation for the effective energy density as a function of 
the total strain. It is as follows: if a strain is inside of the set V, then the energy 
density is asymptotically close to zero as the number of the network nodes goes to 
infinity. If the strain is outside of this set, the energy density is proportional to the 
square of the distance between the strain tensor and V. 

The paper is organized as follows. Section 2 contains formulation of the problem 
and derivation of compatibility conditions for arbitrary deformations. Small deforma- 
tion compatibility conditions are also obtained here. In Section 3 we describe several 
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explicit constructions of still states for arbitrary deformations, as well as a set of 
eigenstrains of all still states. The remainder of the paper deals with the small de- 
formation case. In Section 4.1 we formulate the problem. A new formula for average 
strain is obtained in Section 4.2. Unlike existing formulas that relate average strain 
to nodal displacements, our formula gives an expression in terms of rod elongations. 
This allows to relate components of strain to the average elongations along lattice 
directions. Such formulas seem to be of independent interest. In section 5 we pro- 
vide a complete analysis of the compatibility conditions. In Section 6, we prove a 
characterization theorem for V. This is done by constructing another family of still 
states called stripes. The eigenstrains of the stripes fill T) densely, in the sense that 
every strain within this set can be approximated by an average strain of some stripe, 
with an error of order iV~^/^ where N is the number of nodes. Section 7 contains a 
discussion of the effective energy density. Conclusions are provided in Section 8. 

2. The problem. 

2.1. Energy of one bistable rod. Consider a bistable elastic rod with the 
energy Wr{x) that possesses two equal minima Wr{l) = Wr{al) = 0. Let us also 
assume that Wr{x) is convex outside of the interval [/, /(I + s)], and its second 
derivative is positive in a proximity of each minima x = I and x = l{l + s). For 
definiteness, we may assume that the dependence of Wr{x) on the length x > is 
piece- wise quadratic (W^) or polynomial (WP) 

W,^{x) = min |i(.T - 0', lc{x - 1{1 + sjf 

WP{x)^{x^ir{x-l{l + syf, (2.1) 

where I is the length of the rod in the reference configuration, / -I- s is the length in 
the elongated mode, and i < 1 + s < 2. 

The rods are elastic and locally stable in a proximity of the equilibria. The elastic 
force fr in the rod is 

... dW, 

The rod has two equilibrium states x = I and x ~ 1(1 + s) of equal energy 

fr{x)^0, Wr{x)^0 iix^l, l{l + s). 

The magnitude of the force monotonically increases with the elongation I in the prox- 
imity of equilibria. 

There are several equilibrium lengths Xa and xp for every fr in a proximity of 
zero. They are solutions of the equation 

dWr 

dx 



2.2. Triangular network. Consider a triangular periodic network made of the 
links defined above. For each node, there are six rods joining it with six neighbors. 
When the rods are of equal lengths, the network is elastically isotropic, its Poisson 
coefficient is equal to j, [8]. The number of nodes in the period is three times less 
than the number of rods between them. There are three families of the parallel rods 
in the network 
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I l{l + s) 



Fig. 2.1. Energy of a bistable rod 




Fig. 2.2. To compatibility 



When the rods transit to a different state, the network becomes inhomogeneous 
(each rod can have a different length). We assume that the transition is periodic. Each 
period consists of N nodes where N can be arbitrary large. The network's energy 
is the sum of energies stored in all 3iV links. The length of a link can be expressed 
through the position of its ends, that is the nodes. The nodes are determined by N 
pairs of coordinates in a plane, or by 2N parameters. Hence, the links' lengths can 
not be arbitrary (for example, see Figure 1^?^ . We conclude that they are subject to 
N compatibility conditions that we derive now. 

2.3. Compatibility conditions. 

2.3.1. Necessary condition. There are conditions that constrain the lengths 
of the rods that can be joined in a hexagonal lattice. Consider an inner node in the 
lattice. There are six rods that link the node with the neighbors; let us denote the 
lengths of them as ai , . . . ag . Consider also six rods that surround the node forming a 
hexagon around it; let us denote the lengths of them as 6i, . . . 6g. These rods form six 
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triangles joined at the node. The angles at the node will be denoted by 0i, . . . , (/)g . 

The lengths of the listed twelve rods cannot be arbitrary. The sum of angles of 
the six triangles aia2,bi, .. is equal to 27r, 



^</), = 27r. (2.2) 



i=l 

To express this condition in terms of the lengths, we use trigonometry. Cosine of each 
angle (pi is expressed through the lengths of the links as 

cos(0.)= °' t°'^^~^\ . = 1,...,6. (2.3) 

Here, ay = ai. Therefore, the lengths of rods are constrained as follows 

^ , ,,,, 

> arccos = 27r. 2.4) 

The number of these constraints is equal to the number of inner nodes in the structure. 

2.3.2. Linearized compatibility conditions. In a linearized case when the 
lengths of the rods are close to I, constraints (|2.4p are simplified. Assuming that 

= 1,...,6 (2.5) 
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we find 



^ 0j = 27r + — ^ Ka, + ^ Ka, + i - 2 ^ Kb J 
1=1 \i=l 4=1 1=1 / 



Substituting this expression into (|2.2p . we obtain an elegant linearized compatibility 
condition for lattices 

6 6 

1=1 i=l 

It states that the sum of the elongations of the spokes that come out of a node is 
equal to the elongation of the hexagonal rim around this node. 

3. Eigenstrains and still states. 

3.1. Definitions. 
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3.1.1. Still configurations. In this section we consider still configurations. 
These are deformed states with zero force in each rod. Because of the force-elongation 
dependence, this requires the length of each rod to be either I oi 1(1 + s). To simplify 
presentation, in this section we scale / = 1 and denote 1{1 + s) — a. Given a collec- 
tion of such rod lengths, we call a configuration still if it is geometrically compatible, 
that is the ends of the link of the length one and a can meet in the nodes. The 
examples of still configurations are given below in Figure 13.41 This section contains 
somewhat informal, but explicit description of the set of (average) eigenstrains of all 
still states corresponding to arbitrary (not necessarily small), deformations. We recall 
that eigenstrain is a generic name given to various non-elastic strains, such as strains 
due to thermal expansion, phase transformation, initial strains and so on. 

Notice that the problem is reduced to a geometric problem of tiling of the plane 
with triangles of four types with the lengths of the sides being equal to (1,1,1), 
(1,1, a), (l,a, a), and (a, a, a), respectively. 

Any still configuration can be split into parts that are also still configurations. 
Each configuration is characterized by a symmetric second-rank eigenstrain tensor 
which shows the change of the initial shape of a periodic cell. Assume that a unit 
square element of periodicity of an initial (non-transformed) network contains a very 
large number of triangles. After deformation, it changes its shape. Asymptotically, 
this element is close to a parallelogram. The eigenstrain tensor describes the shape 
of it: the lengths of its sides and the angle between them. 

3.1.2. Eigenstrain of still states. Generally, a deformation due to transfor- 
mation is followed by elastic deformation of all links due to compatibility conditions 
(|2.4p . Therefore, eigenstrain of the composite network is not a convex combination of 
eigenstrains of its components that are, in principle, incompatible and lead to addi- 
tional distortions on the boundaries of components. The eigenstrain of a compatible 
configurations of still states is an affine function (the average) of the eigenstrains of 
the phases. It is independent of the structure of the mixture 

-^mixture = X!^'-^*' = = , * = 0, ■ • ■ , ^ (3-1) 

i i 

where is the volume fraction of the phase with eigenstrain Ei. The only role of the 
structure is the selection of compatible states. 
Next, we describe the variety of still states. 

3.2. The homogeneous still states. In the initial network, there are three 
families of codirected rods. The homogeneous still states are the states in which all 
the rods in one family are in either "short" or "long" mode. 

3.2.1. Classification of homogeneous states. There are eight (2^ = 8) ho- 
mogeneous states that consist of periodic arrays of equal triangles and their reverses. 

1. Initial state Sa with the lengths of the links equal to (1,1,1). The elementary 
square of the state (the eigenstrain) is assumed to be an identity matrix, 

= (l ;) (3.2) 

2. The dual state Sp with the lengths of the links equal to (a, a, a). The unit 
square element Ea of the initial state is transformed to 
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Fig. 3.1. The homogeneous still states 

and it is still proportional to an identity matrix. We call Ef^ the eigenstrain 
of the state 5^. 

3. State S.y with the lengths of the links equal to (a, a, 1) that consists of isosceles 
triangles. The square element Ea of the initial state is transformed to 

/tan7 Q\ 

Here, 7 is the angle between the base of the isosceles triangle and its side, 7 
is expressed through a as 

1 

a = 



2 cos 7 



E^ is the eigenstrain of the state 5-y. 

4. Together with Sj, there are two twin states, S!y with the lengths of the links 
equal to (o, l,o) and with the lengths of the Hnks equal to (l,a, a). The 
eigenstrains of these states are obtained from E^ by rotating it to the angle 
of ±60°, respectively. 

E'^ = R^E^R E'^ = REjR'^ (3.5) 

where R is the 60'' rotation matrix. 

5. State Ss with the lengths of the links equal to (1, 1, a) that consists of isosceles 
triangles. The square element Ea of the initial state is transformed to 

( 2 sin S n\ / 2 sin 5 n \ 

V3 ^] = ( V3 ^ ] (3.6) 

Here, 6 is the angle between the base of the isosceles triangle and its side, 6 
is expressed through a as 

a = 2 cos S 
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The angles S and 7 are related as 



cos 7 cos S = - 
4 

6. State Sg is also accornpanicd by the; two twin states, S'g with the lengths 
of the links equal to (a, 1, 1) and Sg with the lengths of the links equal to 
(1, a, 1). The eigenstrains of these states are obtained from by rotating it 
by the angle of 60°, and —60°, respectively. 

E'g = R^EsR E'^ = REsR^ (3.7) 

Remark 3.1. Pairs Sa , Ss are dual: each of them is transferred 

into another by mutually replacing rods of the length 1 and a. 

3.2.2. Vertices of the eigenstrain states set. Consider a set T> of all pos- 
sible eigenstrains. In a three-dimensional space where coordinates are entries of the 
eigenstrain tensor, this set is a bounded region. The eigenstrains of the homogeneous 

states correspond to the vertices of V. Indeed, there the elongation of the each family 
of the rods is extreme, that is the deformation in the directions of the rods' families 
is extreme. The rods cannot be elongated (contracted) in these directions without an 
elastic deformation. Therefore, there exist three trial matrices Etn Etr^ and Etr^ of 
deformation that correspond to impossible change of eigenstrain in the sense that 

E^ + eEtn , 2; = a, . . . , /3, i = 1, 2, 3 

does not correspond to any still state, for any e > 0. 

Consider, for example, state ^c,. Set of impossible deformations Etr contains a 
cone of matrices of the form 

3 

Y^diUtJ (3.8) 

i=l 

here d < is arbitrary. These deformations require contraction of one of the rods. 

In the state Sa such contraction is impossible. The impossible deformations form a 
cone in the eigenstrain space. Similar cones can be established for all homogeneous 
states because the structure cannot be extended (compressed) in the direction of long 
(short) rods. Below we describe the two- and three-phase composites of still states. 
The eigenstrain of any still state composite is a convex combination of eigenstrains 
of its phases, therefore, the two-phase composites correspond to edges of T) and two- 
phase composites - to its phases. Set I? is a polyhedron. 

We also describe a two-dimensional set TLl of eigenvalues (A^, A^) of eigenstrains. 
This set is also bounded and symmetric (includes pairs (A^,Ab) and (Ab,A^)). 
The eigenstrains of the homogeneous states correspond to vertices oi Hl, and the 
differently oriented homogeneous still states S^, S!y and 5*" (and Ss, S'g and Sg) 
correspond to the same pairs of eigenvalues ~ vertices of • However, the eigenvalues 
of a two-phase composite are convex combinations of eigenvalues of components only if 
their vertices are coaxial tensors (the eigen-axes are codirected). Otherwise, they are 
located along a curve that join the vertices. Therefore, Hl is a symmetric curvilinear 
hexagon. 
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(a) an isotropic and (b) two-phase twin (c) two anisotropic phases 
anisotropic phases 



Fig. 3.2. Examples of laminates 

3.3. Laminate structures. There are still states composed from several basis 
homogeneous states. We call two still states Ey and Ez are called laminate compatible 
if lengths of their rods along the layers are the same. In terms of eigenstrains, the 
states are laminate compatible if a vector r (the tangent) exists such that 

T^{Ey-Ez)T^O (3.9) 

that is, matrices Ey and Ez have the same tangent component. All pairs of homoge- 
neous states except Sa and Sfj can be laminated together, along the side of a common 
self-strain, if the states are allowed to rotate freely. We also remark that threads in 
laminate-compatible structures may possess kinks, as in Figure 4. The continuity of 
the angles is not required for compatibility of a lattice, in contrast with continuum 
mechanics. 

3.3.1. Structures from an isotropic and anisotropic phases. States Sa 
and Sg are laminate compatible. The eigenstrains of mixture are 

k 

EaS = fiEa + (1 - fJ^)Es, /i = " 

where is the volume fraction of phase Sa, ka is the thickness of the layer, and kg is 
the thickness of the layer of phase Sg ■ The laminate is shown in Figure 13.41 

In addition, the eigenstrains Ea and Efs are isotropic and therefore their eigen- 
vectors are codirected with eigenvalues of any phase they are mixed with. A structure 
of one of them and another phase leads to average of eigenvalues Ai , A2 of eigenstrain. 
The eigenvalues A^^ , A^^ of Eag are located at the straight line between eigenvalues 
of Ea and Eg. 

Kg = i^K + (1 - m)a?, a^, = ^iXi + (1 - m)a^ 

Similar laminates can be built for £Jq^, Efjg, E^p, and Eg^. The eigenvalues of these 
structures span the intervals between the eigenvalues of the homogeneous states. 

3.3.2. Laminates from two anisotropic phases Sg and S^. Consider the 
compatibility of two anisotropic states Sg. Each of them consists of "long" and "short" 
rods. Assume for definitencss, that tangent r to laminates corresponds to the line 
of "short" links. This line is an axis of symmetry for the state E.y that consists of 
two families of long and one family of short rods. In orthogonal coordinates r, i^, the 
eigenstrain E^ is diagonal, its eigenvectors are codirected with r and v. 

9 



The state Eg (that consists of two famihes of short and one family of long rods) 
can be made compatible with if it is rotated so that its rr line consists of short 
rods. In other words, tensor of eigenstrain Eg must be rotated to a proper angle so 
that its tangent entry becomes compatible with E^ (|3.4p . [.E^]^^ = 1, see compatibility 
conditions p.9p . The rotated matrix is computed from the additional conditions that 
trace and determinant are independent of rotation. The rotated eigenstrain becomes 

;^sin(6) + 2cos(6) - 1 p 
P 1 



-^<5 (rotated) = ' ^ 



where = ■|(2cos(6) — 1)(2 sin(6) VS — 3). Eigenstrain Es^ of the mixture is still an 
average of those in the phases according to (|3.ip . 



(3.10) 



A^-E .5 (rotated) + (1 ^ Ai)-^7 - 

(^sin(6) + 2cos(6)-l)+(l-p)i^ pp^ 



pp 1 



where p is the volume fraction, as before. 

Eigenvalues A^^ of a composite Eg^ are parametrized by p € [0, 1]. They depend 
on p nonlinearly. In the eigenvalue plane, they are located at an arch that spans Tg 
and A^, not at an interval between them. 

3.3.3. Two-phase twins. Besides these structures, there are two-phase twins: 
the laminates of differently oriented phase Eg or E^. These structures preserve the 
trace of Eg and E-y, respectively, and their eigenstrain is as follows 

_ sin(6) + 2 cos(6) - 1 (1 - 2p)p 



V (1 - 2p)p 

where p = k/N, k e [—N,N]. One could check that the eigenvalues of these twins 



are never isotropic. Figure 3.2(b) 



Together, all listed structures form a region ft in the plane of eigenvalues of E. 
The boundary has three components corresponding to structures Sag (an interval), 
Sg-y (an arch), and 8^/3 (an interval) and four vertices, see Figure [3T3l 



3.3.4. Three-phase laminates. There exist three-phase still structures with 
eigenstrains that densely cover the domain H. bounded by the eigenstrains of two- 
phase structures. Curvilinear triangular domain a, (5, 7 is covered with laminates 
from three phases, obtained by laminating the isotropic phase Sa with the laminate 
Sg.y. The tangents to laminates correspond to the sides of unit lengths in all three 
components. The eigenstrain Eagj is 

Eag^ = lyEa + A^-^S (rotated) + (1 " P - J^)-E^ 

where 

a = k = 0, . . . ,N, V = —, p = 0, . . . ,N - k, 

The corresponding eigenvalue pairs cover the part of that lies outside the rectangle, 
three of whose vertices are a, ry, beta in Figure [ 
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Fig. 3.3. Region of eigenvalues of still states 



3.4. Hexagon-triangles-strips. The central triangle a,ri,(3 can be covered by 
mosaic assembly from elements of S'q, Sfj and differently oriented elements Ss, S'g S'^ 
shown in Figure [3^ Notice that Sa and are incompatible, and the assembly uses 
anisotropic elements Ss between them. The assembly is as follows. 

1. First, three parallelograms from Ss, S'g and Sg are formed of the sizes fc, m, 
fc,n2, and k,n^, respectively. The links of the unit length form the base of 
each parallelogram and the links of the length a form the sides (Figure 13.41 



2. Second, these parallelograms are joined by the edges with the equilateral 
triangle with the side ka in between, forming the structure shaped as Y 



3. Third, the obtained figure is copied and translated. The copies are joined by 
edges, leaving empty triangles and hexagons with the parallel sides of ni, n2 
and 713 unit-length elements. 

4. Finally, the triangular parts are filled with Ss and the hexagon parts with 
Sa- The obtained periodic assembly covers the whole plane (Figure [nUl (c)). 

The obtained periodic assembly covers the whole plane. 
The periodicity element consists of 

- three parallelograms from Ss, S'g Sg respectively, of 2ni fc, 2n2 fc, and 2n3 fc 
triangles each, 

- two triangles composed of S^ of k{k + 1) triangles each. 

- one hexagon containing Sa of 2(nin2 + n2n3 + nsUi) triangles. 
Thus the element consists of 



pairs of triangles. 

The eigenstrain £^hex of the assembly is a sum of eigenstrains of its pieces. 



(a)). 



(Figure EH (b)). 



N = k{k + 1 + rii + n2 + na) + nin2 + n2n-3 + "-3"-i 




(3.11) 



^^^^ 





(a) Step 1 



(b) Step 2 




(c) Step 3, the structure 
Fig. 3.4. Hexagon-triangle strips and the assembly 



where 



D = k{k+ l)Ea + ni kEs + nz kE'^ + kE'^ + (mnz + + n3ni)Ep (3.12) 



E 



hex 



k{k + l)Ea + ni kEs + n2 kE'g + kEg + {nin2 + n2nz + n3ni)Eij 
k{k + 1 + ni + n2 + na) + nin2 + n2n3 + nsni 



3.4.1. Asymptotics. In particular, i^hex is isotropic if ni = n2 — ns = n. 
Then 



-®hex — ^is-^j 



/,;(A- + Ij + /v Tr [E,) + 3n'^a 
k{k + 1) + 3n + 3n2 

12 



(3.13) 



When n — > oo and ^ varies in [0, oo], the isotropic deformation varies from one to 
a. 

The described structures degenerate into two-phase structures E5/3 and EaS- If 
ni ^ n2, ns, k, then 

+ (n2 + n3)£;/3 
-C'hcx ^ — r~; ; ; = E5i3 

If 712 <C ni, fc and 77.3 ^ 77i, /c, we have 

fc(fc + + 771 kEs _ 

fc(fc+l)+77l/fc ~ 

3.5. Covering of the virhole domain fl. Together, the hexagonal and three- 
phase-laminate structures densely fill the whole curved hexagon of eigenvalues of still 
states. The two types of covers are met on the line a, 7 when fc^ — > 1 in hexagons 
and fi —^ in three-phase laminates. AqA^A^A^. The covering is non-unique: For 
instance, one can replace Sa with S'^g and Ss with S*^ everywhere (complementary or 
dual replacement) and obtain another complete coverage. The density of the coverage 
is easy to estimate through the explicitly given parameters k,p. The distance from 
an eigenstrain of a discrete N x N periodic structure to the arbitrary point of the 
domain is of the order of A^^^, as follows from the above explicit formulas. 

4. Small deformations. 

4.1. Linearized elongations. In the remainder of this paper we consider the 
case of small deformations which means that rod elongations are linearized near the 
reference configuration. The linearized relative elongation Ky of the rod (i,j) can be 
written as 

(4.1) 



where q^^j is the unit direction from node i to node j, Ui and Uj are the displacements 
of nodes i and j, respectively. 

Consider a triangular periodic network of N nodes, where N is finite, even though 
it can be arbitrary large. We also suppose that the nodes are arranged into a hexagonal 
shape containing 77 nodes along a side (see Figure 5.1 for an example). 

Factoring out P in (|2.ip and using nodal displacements instead of the rod length, 
we write the energy of a rod as 

CP 

w{K^j) = — min{K^^-, {Kij - sf}, (4.2) 

where, s > is a dimensionless parameter characterizing the critical relative elonga- 
tion. The graph of the dependence w{Kij) is shown on figure HTT] (Compare with (|2.ip . 
where x denotes the actual rod length). 

The displacements Uk in the equilibrium state are found by minimizing the total 
energy of the network: 



W{ui, . . .Mjv) = X w{Kij) = 77; 



Ui 
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Fig. 4.1. Energy of a single rod 



4.2. Average strain. An effective homogeneous deformation state of the net- 
work is characterized by an average (hnearized) strain tensor E. In this section we 
obtain a formula for the average strain in terms of elongations . 

Denote the number of nodes in the network by N . Given a set of displacements 
of nodes iti, . . . mat we define a continuum deformation u{x) that coincides with Ui at 
each Xi. This can be done by interpolating. In the present case, the lattice forms a 
triangulation of the physical domain f2. Therefore it is easy to construct a piecewise 
linear interpolant, using finite elements. After this is done we obtain a continuous 
function u{x),x G il such that 

• u{xi) = Ui for i = 1, . . . iV and Xi being the position of node i; 

• u{x) is linear on every elementary triangle of the lattice. 
To u{x) we associate strain tensor e: 



There exist formulas for E in terms of displacements Ui (see [1], However, for 
our purposes it is better to work with the relative elongations Kij because we would 
often prescribe elongations, rather than displacements. To use standard formulas one 
would have to solve the system (|5.ip (see the next section) that relates elongations 
and displacements. We choose a more direct route that involves relating three average 
elongations along the lattice directions with the three independent components of the 
average strain. 

Rewrite (|4.3|) as 



2 



Average strain E is defined by the formula 




(4.3) 



where is the domain of u and jilj is 
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where T is the set of elementary triangles of the network, |T| is the number of ele- 
mentary triangles in the network, and £a is the strain tensor on elementary triangle 
A: 



£A = 



VttA + (Vma)^ _ 1 



A| Jn 



£{x)dx 



Since the average strain tensor £^ is a 2 x 2 symmetric matrix, it has three inde- 
pendent components: 



a b 
b c 



Consider three lattice direction vectors q^, and (see figure 4.2(a) I 



1 


'2' 


1 


r 


1 


■ -1" 


^^^2 





,92=2 


V3 


' 93=2 


V3 





M2 9l 

(a) Three lattice directions (b) Definition of vector fc^ 

Fig. 4.2. 

We can find the components of E using the values q^ ■ (Eq^), r ~ 1, 2, 3. Hence, 
consider 



Using linearity of u on triangle A for edge of the triangle we find 



(4.4) 



I 



Given a triangle A with vertices h,i,j (see figure 4.2(b)) define vector feA G 



as: 



fcA = 



Khj 
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Using this notation we rewrite (j4.4p as a linear system for finding the components 
o, b,c of E given elongations Kij: 



Q 



a 
b 
c 



(4.5) 



where the matrix Q is given by 

Q = 



'ill 
4,1 



2'?l,l9l,2 
2^2, 192, 2 
2^3, 193,2 



Aer 



9?,2 





' 1 





o" 




1 


V3 


3 




4 


2 


4 




1 


V3 


3 




4 


2 


4 



(4.6) 



and function m 



^ by 
m(a;) — Xi 



"f 0' 




'0 


1" 




'0 0" 





+ X2 


1 





+ X3 


f 



Function m is simply a linear f-f mapping from the space of three dimensional vectors 
to the space of 2 x 2 symmetric matrices. 
Solving (|4.5p we find average strain E: 



(4.7) 



where 



1 

3„ 



3 

1 ■^2 



Further, we will rewrite the sum over triangles in (|4.7p as the sum over edges (this 
is just the change in indexing). Denote all edges of the array by £, all boundary edges 
of the array by £b and all non-boundary edges of the array by Sj {Ei — £\£b)- Then 

^ ^ — 2 ^ ^ fee — ^ ^ fee, 

AgT e6£ eGffi 

where vector fee G M'^ has the components defined by 

, - _ / Ke, if edge e is parallel to _ i o q 

^"^^-^"lo, else ' ^-^'^'^ 

Next, equation (|4.7p can be rewritten as 



E 



\T\ 



Denote the average elongation by fc: 

fe= ^^fce. 



(4.8) 
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Then we rewrite (j4.7p and obtain the desired formula for the average strain in terms 
of average efongations: 

£; = 2^ m(O-ife) - ^ m (q-' ^ fe J . (4.9) 

This formula is used extensively in the remainder of the paper. The first term in the 
right hand side contains contributions from all interior edges of the lattice, while the 
second term contains the contributions from the boundary edges only. As the number 
of nodes increases, the second term becomes a small perturbation of the first. This 
becomes clear from the following estimate 

M M^"' ^ H < V2^||Q-i|bmax||fc,|b. (4.10) 

Here || • \\f denotes the Frobcnius norm of a matrix, and || • ||2 is the Euclidean norm. 
Note that the ratio 

M 

of the number of the boundary edges to the number of all elementary triangles goes 
to zero as the number of nodes approaches infinity. 

5. Compatibility conditions. In the next two sections, we study the linearized 
version of the problem of finding still states. The linearized formulation is as follows. 
Equations (j4.1|) relating node displacements with link elongations can be concisely 
written as 

RU^K, (5.1) 

where 17 is a tuple of displacements Ui, say (mi, . . . itjv), and /« is a tuple of Kij. The 
number of unknowns in (|5.ip is 2N (since each Ui has two independent components). 
The number of equations is larger, approximately 3N (see Lemma 5.1 below for the 
exact count). Therefore, we have an overdetermined system (the number of unknowns 
is less than the number of equations). As a consequence, (15. ip may be unsolvable for 
some K. The values of k for which (|5.ip is solvable are called admissible. 

5.1. Compatibility conditions characterize range of R. In this section we 
show that compatibility conditions (|2.6p ar necessary and sufficient for solvability of 
(|5.ip . Necessity is clear from the derivation of (|2.6p . so we focus on sufficiency. 

Suppose that we have the hexagonal array with a side of n points (see figure [5TT|) . 

For every hexagon consisting of 6 vertices we formulate the condition (|2.6p : the 
sum of elongations of the outer edges equals to the sum of elongations of the inner 
edges. For example, for figure [01 the equation will be 

Kl2 + K23 + ••• + K56 + ^61 = Kji + Kj2 + ■■■ + Kj^ + 

The collection of all such equations will be called hexagonal equations. 

Next, we need a preliminary result giving the exact count of the number of nodes, 
edges and hexagonal equations. The following lemma is used later to carry out proofs 
by induction. 

Lemma 5.1. The total number of nodes N, the number of edges E, and the 
number of hexagonal equations M depend on n as follows. 

iV(n) = Sn^ - 3n+ 1, £'(n) On^ - 15n + 6, M(n) = - 9n + 7. 
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Fig. 5.1. Example of the network with N = 19, n = 3. A hexagon that forms a hexagonal 
equations shown in bold 



Proof. We will prove the lemma using mathematical induction. 
For n = 2 we have a hexagon that has 7 vertices, 12 edges, and 1 hexagonal 
equation. Indeed, 

iV(2) = 3-22-3-2 + 1 = 7, £;(2) = 9-22-15-2 + 6 = 12, M(2) 3-22-9-2 + 7 = 1 

Thus, the formulas are true for n = 2. 

Suppose that given formulas are true for n = k. Let's prove that they are true 
for n — k + 1. 

When we increase rt = fc by 1 we "wrap" the array by a layer of points. Therefore, 
we add 6k vertices, we allow 6(fc — 1) new hexagonal equations, and add 6(3fc — 1) 
edges. Thus, the number of vertices N* for the size n = fc + 1 is: 

N* = N{k) + 6fc = 3fc2 - 3fc + 1 + 6fc = 3(fc + 1 - 1)2 + 3fc + 1 
= 3(fc + 1)2 - 6(fc + 1) + 3 + 3/i: + 1 = 3(/c + 1)^ - 3(fc + 1) + 1, 

the number of equations M* is 

M* = M{k) + 6(fc - 1) = 3fc2 - 9fc + 7 + 6(fc - 1) = 3(fc + 1 - 1)^ - 3fc + 1 
= 3(fc + 1)2 - 6(fc + 1) + 3 - 3A: + 1 = 3(/c + 1)^ - 9(fc + 1) + 7, 

and the number of edges E* is 

E* = E{k) + 6(3fc - 1) = 9fc2 - 15fc + 6 + 18A; - 6 = 9(/s + 1 - 1)^ + 3fc 
= 9(/j + 1)2 - 18(fc + l) + 9 + 3k = 9{k + 1)2 - I5{k + 1) + 6, 

Since N* = N{k + 1), M* = M{k + 1), and E* = E{k + 1), the formulas are true. □ 
Observe that 

2N-3 = E- M. 

We know that for the mapping R : R2JV _^ have dim(keri?) — 3. This 

follows from the well known results on graph rigidity, in particular theorems on the 
first-order rigidity of triangulations (see e.g. [3] and references therein). Thus, 

dim(im R) ^ 2N - 3 
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The hexagonal equations can be concisely written in the form 



Zk ^ 0, 



(5.2) 



where Z is a matrix, Z : — > ; i-th row of this matrix corresponds to i-th 
equation. 

Lemma 5.2. Matrix Z has full rank M. 

Proof. We will restate the lemma in the form: all hexagonal equations are linearly 
independent. 

We will prove the lemma in its new form using mathematical induction on n. For 
n = 2 we have M = 1. One equation is linearly independent. 

Suppose, the lemma is true for n = k. Let's prove it for n — k + 1. 

Increasing n = fc by 1 we wrap a hexagon of side size k with a layer of points. 
Doing so we add some number of hexagonal equations to those that we had in the 
original hexagon. The hexagonal equations of the original hexagon are linearly inde- 
pendent. Each new hexagonal equation contains an edge that is present in no other 
equation. Thus, all equations are linearly independent. □ 

Since matrix Z has rank Af, we find dim(kerZ) = E ~ M. Thus, 



Necessity of hexagonal equations implies that each admissible k satisfies them. 
This fact can be stated as follows: for any k e imi? we have Zn = 0. This in turn 
implies 



Theorem 5.3. For any vector k G satisfying Zk = there is a unique (up 
to translation and rotation) vector u £ satisfying Ru = k. 
Proof. The conditions (|5.3p and (|5.4p imply 



Thus, for any k e R^ we can find a vector u G R^^. If we neglect rigid rotations 
and translations of the whole lattice (they form the null space of R by the well known 
results on rigidity of triangulations, ([3J and references therein), then this vector is 
unique. □ 

6. Still states and small deformation eigenstrains. In the small deforma- 
tion case, a still state is a collection of nodal displacements Ui) satisfying equations 
(j5.ip with the right hand side k of special form. The components take only 

two values: and s/2, where s/2 is the critical elongation from (14. 2|) . If K,y = we 
call the corresponding edge (i,j) short, otherwise the edge is called long. To construct 
a still state, one could choose a length for each edge (long or short). The resulting 
configuration is accepted if the resulting triangles form a tessellation. Otherwise the 
configuration is rejected. In the case of small deformations, an admissible elongation 
vector K, should lie in the range of the matrix R from (j5.1[) . Equivalently, such k must 
be a solution of the hexagonal equations (|5.2p . 

Denote the set of all still states hy U. To characterize U, one needs to solve the 
following problem. 

Given k with e {0, s}, find U realizing that k, that is, find (/such 
that RU = K. 



dim(im R) = dim(kcr Z) 



(5.3) 



im i? C ker Z. 



(5.4) 



im R = ker Z. 
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This is possible only if k satisfies hexagonal equations from Section [51 Complete 
characterization of U seems very difficult, and is not addressed in this paper. The 
reason for the difficulty is the lack of a convenient structure. It is easily checked that 
lA is not a linear space, or even a convex set. Therefore, a description of this set 
cannot be obtained by linear algebra methods. 

6.1. An algorithm of adding still states. In this section we characterize the 
set of strains that can be well approximated by the average strains of still states. An 
average strain of a still state can be conveniently described by a triple of concentrations 
of long edges parallel to the lattice directions. In this section we propose an explicit 
construction that furnishes a large number of still states. Concentrations of these 
states densely fill a certain unit cube in the concentration space. Our construction 
is based on two observations. First, one needs a binary operation that produces new 
still states by combining two already known still states. Second, one needs a simple 
"building block", that is, a still state that could be easily combined with its slightly 
modified (e.g. translated and rotated) replicas. In view of the linear structure of imi?, 
the convenient operation is summation. However, summation of two still states does 
not always produce a still state. An additional condition that ensures that the sum 
of still states remains a still state is non-overlapping of the long edges (see Lemma 
16. II below): if a particular edge is long in the first state, then it should be short in the 
second state, and visa versa. 

From this we deduce that the building block should have low concentrations of 
long edges, and the placement of these edges should be localized as highly as possible. 
If both of these requirements are satisfied, one can add together shifted copies of the 
building block to generate new still states with higher concentrations. 

Lemma 6.1. The sum Ki -\-K2 of two still states Ki and K2 is a still state provided 
Kl ■ K2 = 0. 

Proof. Since Ki and K2 are still states, their components satisfy 

(/«i)y e {0,s}, {K.2)ij e {0,s} 

for all connected nodes i and j. Moreover, the condition of orthogonality Ki ■ K2 ~ 
says that for each certain pair (i,j) the values (/«i)ij and {1^2)1] cannot equal s 
simultaneously. Thus, 

(/«l).y + {K2)ij G {0,s} 

Also, since Ki and K2 are still states, we have 

Zki = 0, Zk2 = 0, 
where Z is the matrix of the hexagonal system (|5.2p . By linearity Z, 

Z{ki + K2) = Zki + Zk2 = 
Since the necessary conditions hold, the sum Ki + K2 is a still state. 

□ 

Notice that Lemma [6 . II states that sum 

Kl + K2 + ■ ■ ■ + K,n 

is a still state, provided Ki, i = 1, . . . , m have non-overlapping long edges. 
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For a still state in (|4.10p we have Hfcelb < s. Thus 



<V2\^\\Q-%s 



(6.1) 



Theorem 6.2. For any ai G [0, 1], i = 1, 2, 3, and any n > 3 there exists a still 
state K* such that its concentrations a*, i — 1,2,3, satisfy 



\ai-a*\ < ^ 



Proof. To prove the theorem we will use special still states called stripes. These 
states have = everywhere except for one stripe that has non-zero elements in 
one direction (figure [Q]) . We can see that these stripes are still states because their 
components are either or s and they satisfy the hexagonal equations: each hexagon 
has either all zero elongations or one non-zero inner elongation and one non-zero outer 
elongation. 



• • • 



• • • • 

• • \ • • 

• • • 



(a) A horizontal stripe, group 1 



(b) Rotated stripe, group 2 



• • 



• • 



• • 




(c) Rotated stripe, group 3 



(d) Compound still state 



Fi g. 6.1. \6.1(a)\\6.1(b)\ and \6.1(b)\ are "stripe" still states. \ 6.1(d)\ is a still state composed of 

\6.1(a)\ \ 6.1(b)\ and \6.1(c)\ Solid line corresponds to k = s, dotted line is 0. 
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Observe, that for two different stripes Ki and K2 we have Ki ■ K2 = 0. In other 
words, the fong edges of /«i do not overlap the long edges of and thus Ki + K2 is 



a still state (figure 6.1(d) ) 



We have 3 different groups of stripes, Qi, Q2, and Q^, one for each direction: one 
group of horizontal stripes and two groups of rotated stripes. The number of stripes 
in each group is 2(n — 1). 

Next, define the vectors dr, r G {1, 2, 3} as follows. The dimension of dr is E (the 
number of edges), and the components are equal to one if the corresponding edge is 
parallel to q^. (a lattice direction vector), and equal to zero otherwise. It is easy to see 
that if K G Qr for some r G {1, 2, 3} then k ■ dr < 2n, r = 1, 2, 3, since each "stripe" 
has no more than 2n long edges. Also, since each stripe has long edges only in r-th 
direction, we get 



This implies 



K ■ dr = — ar, K ■ dt = Q for t ^ r. 
3 



6n 6n 1 n o 

ar < — = 7-^5 7- < — for n > 3. 

E 9^2 - 1571 + 6 71 



First, we explain the main idea of the proof. Consider i-th group (stripes of one 
certain direction). If we add all stripes of this group then we'll get ai — 1. For the 
still state k = Q we have = 0. Therefore, adding stripes one by one we increase 
from to 1 with steps less than n^^ . Hence, for any x £ [0, 1] there exists a step when 
|ai — a;| < n^^. Finally, since different directions are independent from each other, 
we can run this algorithm for i — 1,2,3 and at the end add all still states together to 
get what we need. Next, we present the technical details. 

For any subset S of set Qi we will define functions k{S) and aj(S) by 

Since k{Qi) — di, we have ai{Gi) = 1 and aj{Gi) = for j 7^ i. Also, we define 
k{9) = and aj (0) = 0, j = 1,2,3. 

From the corollary of lemma lOl it follows that for any S C Gi the value of function 
k{S) is a still state with concentrations C(j{S). 

Suppose we are given values Ui £ [0, 1], i = 1, 2, 3. Let's define sets Hi, i — 1,2, 3, 

as 

n^^{Sc g., : a,{S) < a,} 
Having defined Hi we define values S* , i = 1, 2, 3, by 

S* — argmaxdi(5) 



Since the sets Gi are finite, the sets Hi are finite and the maximum exists. Therefore, 
the definition of S* is consistent. 

Now if for some i € {1,2,3} we have g,\S* ^ then S* = g^ and a{S*) = 1. 
From the definition of S* we have 



ai{S*) < ai 
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Since < 1, we get ai{S*) — a^. 

If for some i G {1,2,3} we have Gi\S* ^ then there exists k £ Gi\S* 
From the definition of function a we get 

aiS* U {k}) = aiS*) + a{{k}) > d(5*) 

And, since S* has maximal value in Hi, we have S* U {k} ^ Hi- Hence, 

ai < ai{S* U {k}) 



Now using 



and 



a^(5*)<a^<a,(5*U{fc}) 



a,{S* U {k}) - a,{S*) = a,{{k}) < - 

n 



we get 



< a* - d,(5*) < a,{S: U {k}) - < - 

n 

Finally, for the still state k* defined by 



we have 



\a* -a^\<^ 

□ 

Let us define a set V: 

V = {sm{Q-^x) : a; e M^0 < a;^ < l,r = 1,2,3}, (6.2) 

where is the inverse of the matrix given by Eq. (j4.6p . Set I? is a convex polygon 
in the space of 2 x 2 symmetric matrices, since it is an image of a cube under a 
non-singular linear transformation. The components of the elements of T) form a set 
shown on figure 16.21 

Theorem 6.3. For any E £ V there exists a still state of the hexagonal array 
with a side of n > 3 points such that average strain E* of this still state satisfies 



£^I|f<s||Q"'||2 



n- 1 



Proof. Since E £1), there exists a E , < ar < 1, such that 

E = siaiQ^^cy:) 



Using theorem 16.21 we can claim that for n> 3 there is a still state with the concen- 
tration values a* such that 

IK-a||2< — 
n 
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(a) The cube [0, if 



(b) Q-Mo,i]' 



(c) Eigenvalues of E 



Fig. 6.2. The cube and its image under Q ^ 



\E* -E\ 



Hence, using (14. 9|) we find 

1^1 ni(Q-ia*)-pLm(Q-i ^ kl) - s m{Q-'(x) 
2s \£\ 



<V2\\Q-% 



3 iri 



a — sec 



Next, 



2s_\£\_ 



-a — soi 



< s\\a* ~ cx\\2 + s 



V3 
< —s + s 
n 



1 



2^ 
3|T| 



2\£l 
3|T| 

%/3 



" 2 



Further, using 

|r| = 6(n-l)2, |£| = - 15n + 6, |£ij|=6(n-l) 
and estimation (16.11) we find 



\E* ^E\\F<V2\\Q-^\\2l—^ 



3\T\ 



V3\+V2\^\\Q-%s 



V2s\\Q-^\ 



'V3 \/3 



n 3{n — 1) 71 — 1 



<s||Q-i| 



' n — 1 



□ 

Let us now discuss the significance of the Theorem 6.3 for linking the microscopic 
(lattice- level) and macroscopic (continuum) models. First, it is of interest to under- 
stand how geometry of V depends on the microstructure. We observe that V is the 
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image of a unit cube under a linear mapping defined in (j6.2D . Tlie unit cube and 
the mapping m are both independent of the lattice geometry and particulars of the 
microscopic force definition. The two microstructure-dependent quantities in (j6.2p are 
the relative (non-dimensional) critical elongation s and the matrix Q^^. The actual 
critical elongation is si, where I is the equilibrium edge length. Passing to the F- limit 
requires proper scaling of I with N, a typical scaling being / ~ N~^^'^. In the most 
natural setting, si should scale as /. If that is the case, then s — 0{1) as N ^ oo. 
The set T) in that case is independent of N, and s is a non-dimensional constant of 
the problem that determines susceptibility of the material to phase transition. The 
diameter of D increases linearly with s. The matrix determines orientation of I? 
in the strain space and provides the explicit dependence of V on the lattice geometry. 
Indeed, by definition (|4.6p . Q depends only on the products of components of the 
lattice direction vectors q,., r = 1, 2, 3. 

If s = 0(1), then V is the "flat bottom" of the macroscopic energy density 
because, by Theorem 6.3, any strain in this set can be approximated by a still state 
eigenstrain and the error of this approximation vanishes in the limit N oo. Since 
the energy of every still state is zero, and V is fixed as iV oo, the conclusion follows. 

7. Effective energy of the network. In the previous section, we described 
the set V of all strains that corresponds to zero effective energy. Here, we suggest the 
effective energy equation that corresponds to strains outside this set. 

The average elastic energy W of the network without the transition is described 
by Cauchy formula [8] 

Tr(£2)_i(Tr£f , 

where e is the average strain of the network. Assume first that the transition state is 
fixed in one of the still states E Q V. In other words, assume that no Long-to-Short 
and Short-to-Long transition are allowed. Let us denote the total average strain as e. 
This strain consists of an eigenstrain E and elastic strain e: 

e = E + £. 

The network is elastically inhomogeneous. However, we will assume for simplicity 
that the elastic properties stay constant. We neglect the effect of density variation 
in the transition which is reasonable for small elongations s. Under this assumption, 
the average elastic energy of the network is 

W{£) = W{e-E). 

Minimizing the energy over all still state strains E, we conclude that, for each fixed 
n, the energy density J„ can be written as follows 

iJe) = miiiW(e- E) =0{n-^) ifeeV. (7.1) 

Assume now that e ^ T). The density becomes 

JJe) ^ imnWle- E) (7.2) 

Notice that J is a convex linear combination of the energies of homogeneous states 
that correspond to the vertices of T). This follows from convexity of V and W{e). 

25 



Motivated by this, we suggest a representation for the effective energy density J in the 
hmit n — > oo (r-Hmit). By the well known argument [5], using density of piecewise 
linear deformations, it is enough to consider only linear deformations (equivalently, 
constant effective strains). For these strains, 



Here Ei are the eigenstrains corresponding to the corners of T). 

8. Conclusions. We studied the discrete model of solid-solid phase transitions. 
The model consists of a two-dimensional triangular periodic lattice of nonlinear bi- 
stable rods. In order to characterize the effective energy of the lattice, we obtained 
the description of the set of the average eigenstrain tensors corresponding to so called 
still states, e.g. microstructure deformations that carry no forces. We also proposed 
several explicit constructions of such states for small and large deformations. 

The compatibility conditions for rod elongations were derived. In the small de- 
formation case, these conditions are necessary and sufficient for the existence of a 
deformation realizing a given set of elongations. To characterize the states with zero 
effective energy, we proposed a construction of a special class of still states. The 
eigenstrains of these special patterns densely cover the set T) of all eigenstrains corre- 
sponding to still states. We showed that the set V for our network is a parallelepiped 
in 3D space of independent components of the strain tensor. The orientation and side 
lengths of T) were explicitly described in terms of the parameters of the underlying 
microscopic discrete model; the side lengths are proportional to s and orientation is 
determined by the components of the lattice direction vectors via the matrix Q^^ . 
Components of the inverse matrix Q are defined via components of the direction 
vectors according to (|4.6p . 

The above results suggest the following properties of the energy density: 

• if a strain is inside of T), then the energy density is asymptotically close to 
zero as the number of the lattice nodes approaches infinity; 

• If the strain is outside of this set, the energy density is proportional to the 
square of the distance between the strain and T). 
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